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A proposed solution of the anomalous behavior of the electron spectrum near the endpoint of 
J-^ ' tritium /3-decay is offered, ft is based on a new theory of mass in which mass becomes a dynamical 

Q variable, and the electron in the tritium /3-decay has a narrow mass distribution. The predicted 
Kurie plots explain the main feature ("mt < 0") of this anomalous behavior. 
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Tritium beta decay shows anomalous behavior, observed in many experiments since 1991 Though various 

theoretical explanations have been offered flQ-O], there is as yet no consensus on the solution of the puzzle. The 



anomaly consists of two effects, (a) The Kurie plot lies above a certain straight line (the theoretical result for a 
massless neutrino) near the endpoint, and overshoots that endpoint This is the effect that is fit by the infamous 
"negative neutrino mass squared" parameterization, (b) There is a narrow and low "bump" in the Kurie plot starting 
O a few (5-10) eV below the endpoint Q |. 

We show that if the electrons emitted in tritium /3-decay have a narrow mass distribution instead of the sharp mass 
qj mo assumed in the present day standard theory, the endpoint-ovcrshoot effect (a) above, can be explained. This idea 
is based on a theory of mass developed by one of us; a few details will be given below. The electron antineutrino may 
be massless or have real mass m u > 0, though the latter case gives a better fit, see later. The half-width Am e of the 
electron mass distribution must of course be much smaller than the central value mo (~ 1/2 MeV, units: % = c = 1), 
^ ! and the fits below suggest Am e « m v . 

It is easy to see qualitatively why such a mass distribution would imply the endpoint-overshoot effect. At a value 
of electron momentum p far from its endpoint value the distribution would produce a narrow scatter of events with 
electron energies both above and below Eo(p) = \/ p 2 + TOq • But near enough to the endpoint p = p max (Definition: 
Pmax is calculated from energy conservation for antineutrino momentum q = and mass m v > on the standard 
theory, electron mass sharp = mo ) this scatter would become unsymmetrical because values of electron mass m > mo 
would be ruled out by energy conservation. At p m ax the electron count goes to zero on the standard theory. But 
there would still be counts on the mass distribution theory due to values m < mo . Correspondingly, the new Kurie 
plot K(p) would coincide with the old one Kq(j>) for p far from the endpoint, but at p = p m ax , where Ko(p max ) = , 
K(p m ax) would be positive. Thus K(p) would overshoot the endpoint by an amount of order Am. Figure 1 illustrates 
this quantitatively and in detail. 

To avoid some obvious misunderstandings we add a few words here on the theory |l~5] , |l(| underlying this proposed 
solution. Some further details will be given in the concluding remarks. 

The unsharp mass of the electron (the mass distribution) is not due to decay channels — the electron is stable. 
Rather, mass m of any elementary particle becomes a new dynamical variable alongside linear momentum p and 
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energy E. Mass m is conjugate to a new length variable A just as p and E are to position r and time t, and roughly 
obeys the uncertainty relation Am AA > 1/2 in any state. The new dimension A enters this theory as a microscopic 
length which modifies the usual point-particle causality of the A-D theory at small distances ]l5|. It prevents the 
instantaneous action of source point on field point, and can be thought of as giving (massive) elementary particles a 
finite size ~ A JL7| |. The uncertainties Am and AA depend on the electron's state, here a free electron state. Thus 
Am here is expected to be different from the Am for, say, an electron bound in an atom or molecule. 

This theory is invariant under the Poincare group in particular, so that momentum and energy are conserved. 
Further, it implies the usual relation E = \/p 2 + m 2 for free particles. 

Traditionally in physics every particle had a definite (sharp) mass in the absence of interactions leading to decays. 
But both theory and experiment in the last fifteen years suggest that we may be in a transitional stage in our 
understanding of elementary particle mass. The concept of "particle mixing" has entered particle physics. What 
this implies is that there exist particles, some observable, which have no definite mass. (We mean here: apart from 
the mass widths implied by possible decay channels.) Examples are the three kinds of neutrinos, K° and K°, the 
"primed" quarks d! , s', and b', which directly enter the weak interaction Lagrangian, etc. These particles "oscillate": 
their states are superpositions of a few other particle states with sharp mass, the "mass (eigen) states" . In other words, 
the oscillating particles have (simple) mass distributions. 

However, there are inconsistencies and lacunae in this transitional theory of mass. For example: (i) Among 
elementary particles observable as free particles there are the heavy boson Z° and the photon 7, which are mass 
eigenstates, and then there are the neutrinos v e , and v T , which are mixtures, (ii) Further, though the term "mass 
(eigen)state" has entered the literature (jl8|, the corresponding "mass operator" for the general particle state has never 
been defined, (iii) There is no theory for the mixing angles, nor for the discrete and finite mass spectra we see in 
nature, (iv) Why are massless neutrinos necessarily only left-handed? 

Suppose then that the electron has a narrow mass distribution P(m). Whereas the usual free (Interaction Picture) 
quantum electron field appearing in the 5-operator is |ll| 



i/jq(x) = (2ir) 3/2 J d 3 p (mo/ E Q ) 1/2 e lp ' x u(p,m ) ffl(p) + antiparticle part, 
the corresponding 4-D field with the mass distribution amplitude A(m) is 

tjj(x) = (27r)~ 3 / 2 / / dm d 3 p (m/E) 1 ^ 2 A{m) e lp ' x u(p,m) a(p, m) + antiparticle part, (1) 



where 



E(p,m) = \f p 2 + m 2 , (ij-p + m)u(p,m) = 0, [a(p,m), (p,m)] + = S(p — p')5(m — m'), 



and where P(m) = |A(m)| 2 and / °° dm\A(m)\ 2 = 1. Note that both ipo(x) and 4>(x) have dimension L -3 / 2 (L = 
length) as they must to give a correctly normalized interaction Lagrangian. 

Assume now that the electron field, Eq. (|l]), appears in the S'-operator. Then on taking the S'-matrix element 
corresponding to the decay 3 H — > 3 He + + e + P, the electron differential spectrum J\f(p, m) (where, by definition, 
M(p, m) = number of electrons per unit time with momentum p and mass m in the intervals (p,p+dp) and (m, m+dm), 
respectively) becomes [ p0| 

Af(p,m) = N(p,m)P(m), N(p,m) oc p 2 qE u . (2) 

Here, p = \p\ is electron momentum, q — \q\ is antineutrino momentum, E v = \J q 2 + m 2 is antineutrino energy with 
a (real) antineutrino mass m v > 0, and m is electron mass. 
Energy conservation reads 

Q = M{ 3 H) - M{ 3 He+) = E e + E v , (3) 

where (as usual) we neglect the 3 He + recoil energy and for simplicity treat the case of a single 3 He + final state. In 
the NR (nonrelativistic) case valid for tritium /3-decay we can expand the electron energy as 

E e = [p 2 + (m +2;) 2 ] 1 / 2 w mo + T + x, x = m — m , T = p 2 /2m , (4) 

correct to terms linear in x (assume P(m) « unless \x\ « mo). If Eq. (^) is used in Eq. (j|), NR energy conservation 
reads 
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T = T + x + E v , T = Q - mo. (5) 

We remark here that To is usually called Eq (w 18,570 eV) in the experimental literature. Note that T is not the 
kinetic energy of the electron of momentum p and mass m because mo, not m, appears in the denominator. But it is 
convenient to use this T because then the full m-dependence of E e to 0{x) appears in the single term +x in the NR 
case. Further, T is the kinetic energy in the standard theory, and is what is used in the experimental papers (where 
it is usually called E). 

After putting q and E v in terms of p and x, we get the electron differential spectrum in the form 

Affam) cx 2m T(T - T - x) [(T - T - x) 2 - ml} 1/2 P(m). 

It is more convenient to use a spectrum Af'(T,m) for the number of electrons with kinetic energy T between T and 
T + dT, etc., so we multiply the above spectrum by dp/dT — m§j \J2moT . Changing to the convenient energy variable 
y = To — T , we arrive at 

M\y,x) cx m y/2m (T a -y) (y - x ) [ ( y - x f - ml] 1 ' 2 P{x) dydx, (6) 

where we wrote Af'(y,x) = Af'(T,m) and P(x) = P(m) for notational simplicity. Note that since momentum p (or 
equivalently, kinetic energy T = p 2 /2mo ) is measured in the experiments but not mass m, the differential spectrum 
of interest is 

A/-'(y) = / N'(y,x)dx. (7) 

Integrated spectrum. In the experiments , when a Kurie plot is given, it is the integrated spectrum (differential 
spectrum integrated in T from some T to the endpoint). Hence we consider the integral of Af'(y) from the endpoint 
y = to some y > 0. The only values of y of interest are in the endpoint region y << To, so we can put To — y — * To 
in Eq. @. Then from Eqs. (§) and (@): 

/V—m v ry 
dxP(x) dy'(y' - x)[(y' - x) 2 - ml] 1 ' 2 . (8) 

-nAm J x-\-m v 

The limits were determined from energy conservation, Eq. ([5]), on the assumption that the y'-integral is done first for 
a given value of x. In the x-integral, Am is the half-width of the mass distribution, and n is some small integer (2 or 
3) such that P(x) is considered effectively zero for x < —nAm. 

Remarkably, the ^/'-integral in Eq. @ can be done analytically. Then we define a normalized pi[ integrated Kurie 
plot K (y) by dividing out some factors and taking the cube root of the resulting integral I(y) : 

dxP{x)[{y - x) 2 - ml] z/2 , m v - nAm < y « T , (9a) 

-nAm 

K(y) ee lliy)} 1 / 3 . (9b) 
The integral I(y) must be done numerically, and we chose P(x) to be the gaussian 

p{x) = 7ib cxp (- (10) 

for the numerical integration, in which case we set Am = a, the standard deviation of the gaussian. But see the 
concluding remarks. 

There are two mass parameters available, a and m v . It is convenient to non-dimensionalize Eqs. (^) and ( |l0| ) by 
using one of them. We chose m V) and denoted the dimensionless quantities by carets: 

x = x/m v , y = y/m v , d = a/m y , T = T /m„; P{x)=m v P{x), I(y) = I(y)/ml, K(y) = K(y)/m v . 

Then the quantity plotted in Fig. 1 is K (y) = [/(y)] 1 / 3 , where 

T{y) ee j V dxP(x)[{y - x) 2 - l] 3 / 2 , l-na<y«%, (11) 
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for various values of the ratio a = a/m v . 

Note that in the sharp mass limit P(x) goes to a delta function, I(y) becomes (y 2 — 1) 3 ^ 2 , and 

K{y) -> (y 2 - 1) 1/2 « y - 1/2% for y » 1. 

This is a straight line of slope +1 (in y ) far from the endpoint, but near the endpoint it curves down and vanishes 
at y — 1 below the endpoint, the well-known result for a massive neutrino in the standard theory. (For m v — , the 
corresponding plot of K(y) would be a straight line all the way to y — 0.) 

In the present, mass unsharp theory for the case m u — we must return to K(y), Eq. (|9b|) , with m u set to zero. 
This can be non-dimensionalized by dividing all quantities by suitable powers of a. We have omitted this graph here; 
the Kurie plot lies approximately half way between the plotted curves for a = 1.0 and 1.5 shown in the figure. 




with experimental results), for a = 0.3, 0.7, 1.0, and 1.5 from left to right. Also shown are the Kurie plots of Ko{y) in the sharp 
mass case for massive and massless neutrino. 



Discussion of the results. Fig. 1 shows the non-dimensionalized Kurie plot for the cases a — 0.3, 0.7, 1.0, and 1.5. 
For ease of comparison with experimental results, however, we have plotted — y = (T — To)/m v on the horizontal 
axis. Also shown are the Kurie plots in the sharp mass case: Ko(y) — (y 2 — l) 1 / 2 for m„ > and the straight 
line Ko(y) = —y for m v = (the latter is the dimensional Kurie plot Ko(y) — —y arbitrarily normalized with an 
m v 7^ 0). K(y) stays above the straight line and overshoots the endpoint y = (T = To) by about 2mi,,3m„, or 
4m„ in the cases a — 0.7, 1.0, and 1.5, respectively. If one knew or guessed a value for m u , this would give the 
predicted overshoot quantitatively. The comparison with experiment would be clearer if the data at and just beyond 
the endpoint were cleaner, as previously mentioned The curves for the values a < 0.3 lie below the straight line 
and do not overshoot the point T — To. There is no indication in any of these curves of effect (b), the "bump". The 
plots in the Figure resemble the experimental plots in Refs. [1,2,3,5,8, and 9] in a general way. In particular, the plot 
for it = 0.7 appears to be an especially good fit to the experimental points in Fig. 1 of the recent work [9]. 

Concluding remarks, (a) The gaussian, Eq.(10) was chosen arbitrarily as a typical narrow mass distribution, but 
calculations using other mass distribution functions indicate that the Kurie plots do not depend sensitively on the 
shape of the distribution so long as it is narrow. We should mention that the underlying theory ]l5[ | is capable of 
predicting this mass distribution uniquely. Further, recent work jl6| suggests how the effective 4-D field ip(x), Eq. 
(|l|), arises from the natural 5-1? free quantum field ^(x, A) of this 5-D theory. However, these ideas are not yet in 
final form and anyway would be too lengthy to give here, (b) The possible objection that the uncertainty in the 
electron mass found in particle property tables is too small to explain this 3 H /3-decay anomaly is no valid criticism of 
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this proposal. To repeat: if mass becomes a dynamical variable, its uncertainty depends on the state, and may vary 
wildly with that state. The value 0.3 ppm for the uncertainty quoted in the Tables j2^] comes from the most precise 
atomic measurements. There is no reason for the mass width to be the same as that for a free electron especially 
if this "free" electron is acted on by a self-force. (A self-force was found essential in the classical theory to give a 
nonsingular self-interaction with all the correct properties [ fL5| and in the quantum theory to enable the prediction 
of mass spectra for elementary particles |l6|]). To our knowledge no experiment measuring the mass width of a free 
electron, such as emerges in tritium /3-decay, has ever been performed, though it would be simple in principle, (c) We 
can give an upper bound on a = Am e on the basis of our results as follows. We predict a = am v w 0.7m v (see 
above). If m„ < 1 eV for example, then a < 0.7 eV, or < 1.4 ppm. 
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